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Abstract. We study the formation of neon-enriched donor stars in ultracompact X-ray binaries (P or b < 80 min) 
and show that their progenitors have to be low-mass (0.3 - 0.4) Mq "hybrid" white dwarfs (with CO cores and 
thick helium mantles). Stable mass transfer is possible if in the initial stages of mass exchange mass is lost from 
the system, taking away the specific orbital angular momentum of the accretor ("isotropic re-emission"). The 
excess of neon in the transferred matter is due to chemical fractionation of the white dwarf which has to occur 
prior to the Roche lobe overflow by the donor. The estimated lower limit of the orbital periods of the systems 
with neon-enriched donors is close to 10 min. We show that the X-ray pulsar 4U 1626-67, which likely also has a 
neon-enriched companion, may have been formed via accretion induced collapse of an oxygen-neon white dwarf 
accretor if the donor was a hybrid white dwarf. 
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1. Introduction 



Recently, Juett et al. (2001 ) reported the discovery of 
the excess of neon in four low-mass X-ray binaries: 4U 
0614+091, 2S 0918-549, 4U 1543-624, and 4U 1850-0857. 
It is suggested that the Ne is local to these objects. 4U 
1850-0857, which belongs to a globular cluster, has a mea- 
sured orbital period of 20.6 min. From the optical prop- 
erties of the systems, Juett et al. suggest that all of them 
are ultracompact: P or h < 80 min. A local excess of neon 
is found also in the ultracompact (P or b = 41.4 min) bi- 



model and related problems are considered in Sect. [|. Our 
conclusions follow in Sect. [|. 



2. Evolution of ultracompact binaries 

2.1. Formation of ultracompact binaries 



work of Paczynski (1967) on the 
(~ 18 min) cataclysmic variable AM 



nary pulsar 4U 1626-67 (Angelini et al. 1995; Schulz et al 



2001). Juett et al. propose that in all these binaries the 
donors are low-mass, neon-rich degenerate dwarfs. [Bchulz 



et al. (2001) suggest that the donor in 4U 1626-67 is the 
core of a carbon-oxygen (CO) or oxygen-neon (ONe) white 
dwarf which has previously crystallised. However, the evo- 
lution of neutron star+white dwarf [henceforth, (wd,ns)] 
binaries to the Ne-enrichment stage was not studied. 

Below, we explore the evolution of white dwarfs which 
fill their Roche lobes having neutron star companions. 
We discuss the mechanism of mass exchange and limit- 
ing masses of the donors in these binaries in Sect. |^. The 
process of chemical fractionation that could lead to Ne- 



enrichment is briefly outlined in Sect. 2.3. The observed 
systems are discussed in Sect. 0. Some uncertainties of the 
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After pioneering 
ultrashort-period 
CVn, it is commonly accepted that the donors in semi- 
detached systems with the orbital periods of several tens 
of minutes may be degenerate dwarfs and that the driving 
force of the binary evolution is angular momentum loss 
(AML) due to gravitational wave radiation (GWR). 

The basic features of the scenario for the formation of 
(wd,ns) systems which can evolve into ultracompact X-ray 
binaries may be summarized as follows (see e. g. Tutukov 
|fc Yungelson 1993 : [ben et al. 1995, for details). In the 
Galactic disk stars with masses > 10 M Q that are born 
in binaries evolve into neutron stars, experiencing under- 
way one or two common envelope stages which strongly 
decrease the orbital separation. Next, their companions 
experience dynamically unstable mass loss and become 
white dwarfs. If the resulting systems have orbital sepa- 
ration ~ Rq , AML via GWR may bring the white dwarf 
to Roche-lobe overflow (RLOF) within the lifetime of the 
Galactic disk and under certain conditions (see below) sta- 



2 



Yungelson, Nelemans & van den Heuvel: Formation of neon-enriched donors in X-ray binaries 



ble mass exchange leading to exposure of the Ne-enriched 
core may ensue. 

In globular clusters neutron stars most probably ac- 
quire their companions via exchange interactions with pri- 
mordial binaries (e. g. Rasio et al. 2000). The subsequent 
evolution is similar to that of the disk binaries. 



2.2. Mass transfer in ultracompact binaries: the nature 
of the donors and the limits of their initial mass 

In systems with degenerate donors, at the onset of Roche- 
lobe overflow the mass transfer rate M may be several 
orders of magnitude higher than the Eddington limiting 
accretion rate (M E dd) for a white dwarf or ne utron star 
accretor ( |Vila 197lfc flutukov fc Yungelson 1979|) , [ see Fig. 

Formation of a common envelope engulfing the whole 
system and merger of components or collapse of the neu- 
tron star into a black hole may be avoided if the excess 
of the matter which cannot be accreted leaves the system 
taking away the specific orbital angular momentum of the 
accreting component ["isotropic re-emission" (Soberman 
et al. 1997 )1- This role may be accomplished by radia- 
tively driven outflows (King & Bcgclman 1999; Tauris 
fc Savonijc 1999) ). The equation for the mass loss rate 
by a white dwarf in the presence of AML by GWR and 
"isotropic re-emission" used in this study may be found 
in Yungelson et al. (1996). 

Figure [j] shows this mass loss rate for a (wd,ns) system 
with initial masses of donor and accretor Md = 0.83 M© 
and M a = 1.433 M©, respectively]^]. Since all white dwarfs 
approach zero-temperature radii after cooling for several 
100 Myr ( Panei et al. 2000| ) and since most systems will be 
very much older than this when their orbits have decayed 
sufficiently to start mass transfer, we use in our computa- 
tions the mas s - radius (M — R) relati on for cold spheres 



CO donor+neutron star accretor 




M(M ) 

Fig. 1. The solid line shows the mass loss rate from a white 
dwarf "stabilised" by isotropic re-emission for the case of 
a system which initially contained a 0.83 M© white dwarf 
donor and a 1.433 M© neutron star accretor. Dashed line: 
the upper limit of the mass loss rate from the system by 
the isotropic re-emission (for r ns = 10 km). Dotted line: 
MEdd for a neutron star. Asterisks mark the parameters 
of the system at logf(yr) =3, 4, 5, 6, 7, 8, 9. Systems 
with different initial Md follow curves which are indis- 
tinguishable on the scale of the figure. To the left of the 
beginning of the mass-loss curve (A) the Roche-lobe filling 
white dwarfs are dynamically unstable. Between A and B 
isotropic re-emission is not efficient enough for ejection 
of all matter lost by white dwarf. "Hybrid" white dwarfs 
which are able to transfer matter stably have initial masses 
corresponding to the thick part of the mass-loss line BC. 



de rived by |Zapolsky fc Salpeter (1969[ ) in the form given 



by Ra ppaport et al. (1987 ), assuming equal mass fractions 



dwarfs as donors as these have M ^ 1.1 M© (Gil-Pons & 



of carb on and oxygen. For this M — R relation, if the mass 

retention efficiency of the accretor is close to zero, (wd,ns) 
systems with Md £ 0.83 M© are dynamically unstable. 

All matter in excess of MEdd can be lost from the sys- 
tem only if the liberated accretion energy of the matter 
falling from the Roche lobe radius of the neutron star to 
the neutron star surface is sufficient to expel the matter 
from the Roche-lobe surface around the neutron star, i.e. 
M < M max = M E dd (^L.nsAns), where r ns is the radius 



Garcia-Bcrro 2001). Similarly also massive CO dwarfs are 



of the neutron star ( |King fc Bcgclman 1999 ; Tauris et al.| 
200C ) . The dependence of M max on the mass of the donor 
is also plotted in Fig. [l]. Since the mass-loss line in Fig. [I] 
holds as well for initial Md < 0.83 A/©, within our assump- 
tions "isotropic re-emission" from a neutron star is actu- 
ally possible for donors with initial Md ^ 0.44 M©. Unless 
there is an even more efficient way to stabilise the mass 
transfer than by "re-emission", this excludes ONe white 



excluded as donors. 

The overabundance (relative to solar) of Ne suggests, 
that the progenitors of the donors in the X-ray systems 
under consideration have experienced core helium burn- 
ing including the 14 N(a, 7) 18 F(a, 7) 22 Ne reactions chain. 
Thus the donors cannot be low-mass helium white dwarfs. 

This leaves as the last option the so called "hybrid" 
white dwarfs which have CO cores and thick He mantles. 
They are formed from components of close binaries with 
initial masses in the range from 2.5 to 5 M© which experi- 
ence RLOF prior to core He ignition, become He-burning 
stars and, after completion of core He burning, evolve di- 
r ectly into white dwar f s. Their mass is ~ (0.3 — 0.7) M© 
( |lben fc Tutukov 1985| ; |Han et al. 2000|) . 

In a calculation of the model of the population of 



compact stars in the Galaxy (Nelemans et al. 2001) the 
birthrate of such (wd,ns) systems with M w d < 0.4 M© is 
We neglect the effects of an X-ray irradiation ind uced stel- 2.7 x 10~ 6 yr^ 1 . For comparison, the birthrate of progen- 
itors of X-ray binaries with hydrogen-rich donors in the 
same model is 1.3 x 10 yr -1 . 



lar wind whic h may enhance mass loss by the donor ( Tavani & 
London 1993|). 
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2.3. Neon-enrichment 



CO/He donor+neutron star accretor 



Even after core helium exhaustion the abundance of Ne in 
the cores of hybrid white dwarfs is low, so to observe the 
strong Ne enrichment, the cores have to be crystallised and 
fractionated. These processes take several Gyr, depending 
on the mass of the dwarf, transparency of the outer layers, 
t otal amount of 22 Ne and uncertainties in phase diagrams 
ernanz ct al. 1994 ; Montgomery et al. 1999| ). Figure 
shows that the orbital periods of (wd,ns) pairs enter 
the observed range of ~ lOmin in less than lOMyr. (The 
meaning of the lines and asterisks in Fig. || is the same as 
in Fig. [l]). Hence, the enrichment of the white dwarf core 
by Ne, must have happened before RLOF. 

For an initial metallicity Z=0.02, after completion of 
He-burning, the mass abundance of 22 Ne in the core of 
the star cannot be larger than X22 Nc = 0.02. During crys- 
tallisation of binary mixtures of (C/Ne) and (O/Ne), Ne 
settles in the center and forms a Ne-enriched nucleus which 
contains all 22 Ne formed in the dwarf and has the so called 
azeotropic abundance ^22 Ne of Ne (Iscrn ct al. 1991). 



The mass of this Ne-rich nucleus is 



-A^Ne-r = (X22 Ne /X?2 Ne ) M w d. 



(1) 



If, following Bravo et al. (1992), one assumes X 



2 No 



0.07 - an intermediate value of X^ 2Ne for O/Ne and C/Ne 
mixtures - and takes for M wc j the maximum mass of the 
convective core of a 0.378 Mq "hybrid" dwarf - 0.2 Mq 
(Ibcn & Tutukov 1985), then the mass M^e-r of the Ne- 
enriched core is ~ 0.06 Mq. 

If the white dwarf didn't have enough time to crys- 
tallise, the mass loss uncovers material of the convective 
core of the progenitor of the dwarf. Then the mass fraction 
of Ne may be down to the initial ~ 001 



3. Comparison with the observed Ne-rich X-ray 
binaries 

3.1. 4U 0614+091, 2S 0918-549, 4U 1543-624, and 
4U 1850-0857 

Figure ^ shows the evolution of a semidetached system 
which initially contained a "hybrid" white dwarf and a 
neutro n star, in orbital period - mass loss rate coordinates. 



Takin g] the mass of the Ne-enriched core of 0.06 Mq esti- 
mated in Sect. [2.3] as typical, one obtains a rough lower 
limit of the orbital periods of binaries with the Ne-rich 
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Fig. 2. Orbital period - mass loss rate relation for a sys- 
tem which initially contained a 0.38 Mq "hybrid" white 
dwarf donor and a 1.433 Mq neutron star accretor. The 
meaning of the lines and symbols is as in Fig. [l]. The verti- 
cal dashed lines mark the orbital periods of 4U 1850-0857 
and 4U 1626-67. Within the period range corresponding 
to the shaded area the donors may transfer neon-enriched 
matter. 



2.87 ± 0.16. For oxygen an underabundance is reported: 
0.37 ± 0.06 to 0.52 ± 0.04. Taken at face value, these rela- 
tive abundances translate into a local Ne/O ratio of 0.63 
- 0.77. If true, these values are considerably higher than 
one would expect for the above used abundance of Ne in 
the nucleus of the dwarf ^22 Ne = 0.07. But [sern et al 



(1991 ) notice, that ^22 Nc may be well underestimated by 
a factor ~ 3. If this is true, the Ne/O ratio comes into bet- 
ter agreement with observed values. The mass of the Ne- 
enriched core then becomes slightly higher than 0.02 Mq, 
still comparable with our estimates of the masses of the 
white dwarfs in 4U 1850-0857 and 4U 1626-67. 

The problem of a low predicted Ne/O ratio would also 
hold for initially more massive donors, if stable mass loss 
would be possible for them. In the core of a 1.1 Mq ONe 
dwar f the Ne/O ratio is ~ 0.4 ( Gil-Pons fe Garcia-Bcrro 
2001|) . For a "standard" 0.6 Mq CO white dwarf with an 



initially cquimolar distribution of C/O and traces of Ne, a 
final Nc/O ratio ~ 0.25 is expected in the core ( Segretain 



white pwarf donors of about lOmin. The white dwarf 
enters the range of periods (or masses) for which Ne- 
enrichment is possible in several Myr following the onset 

of RL OF (Fig. B). 

Ou 



ct al. 1994| ) 



All four Ne-enriched X-ray binaries contain, presum- 
ably, weakly magnetised neutron stars. This may b e a re- 
sult of the decay of the magnetic field by accretion ( Faam 



calculations of the mass loss rate also give the fc van de Heuvel 1986 



mass of the donor as a function of the orbital period. 
We thus infer the donor masses for 4U 1850-0857 and 4U 
1626-67 to be 0.027 and 0.01 Mq, respectively (see more 
detailed discussion of the latter system in Sect. 3.2). 

The abundance of Ne relative to solar found in the 
observed Ne-rich X-ray binaries ranges from 1.9 ± 0.3 to 



3.2. 4U 1626-67: an accretion induced collapse? 

4U 1626-67 differs from the binaries discussed above by 
the presence of a 7.7 s X-ray pulsar with a strong magnetic 
field estimated from cyclotron emission: (3.2±0.1)(l + z)x 
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10 12 G, where z is gravitational redshift (Orlandini et al 
1998[)~ Applying a polytropic M — R relation, Vcrbunt 
et al. ( 1990| ) estimated its donor mass as 



4U 1626-67 



0.02 M R . A 



similar estimate of was obtained by Bchulz et al. (2001) 
under the assumption that velocities of the lines in the 
Ne/O complex and their Doppler shifts reflect Keplerian 
motions in the accretion disk. Thus, existing estimates of 
Md are consistent with Ne-enrichment of the donors. A 
very low value of Md in 4U 1626-67 is also favoured by 



o 



the probability of detection considerations ( Verbunt et al 
199(0 

The high magnetic field strength suggests that the neu- 
tron star is young and accretion was negligible. The ab- 
sence of enhanced abundances of O-group elements in the 
spectrum of 4U 1626-67 argu es against the origin of the 



O 



-10 



neutro n star in a SN II event ( Angclini et al. 1995 ). Taam 



0.4 




& van dc Heuvel (1986 ) suggested that this X-ray pulsar 
was formed recently - (2 — 10) x 10 7 yrs ago - by an ac 



M(M a ) 



cretio n indu ced collapse (AIC) of a white dwarf. Vcrbunt 
et al. 199C) rejected the AIC model, based on the ar- 



gument that the time between AIC and resumption of 
mass transfer is ~ 10 8 yr, longer than then assumed e- 
folding decay time (~ 10 7 yr) for the magnetic fields of 
neutron stars. Later, the analysis of magnetic fields of iso- 
lated radio pulsars have shown that magnetic field decay 



scale is probably £ 10 8 yr ( Bhattacharya et al. 199^ ) 



This, along with the discovery of Ne-enrichment in the 
system, suggests that one should reconsider the formation 
of 4U 1626-67 via AIC, especially, since the details of pre- 
collaps e evolution wore never studied before. 



Fig. 3. Solid line - mass exchange rate in a system which 
initially consisted of a 0.38 Mq "hybrid" donor and a 
1.2 Mq ONe white dwarf accretor. The dashed line is the 
upper limit of the "isotropic re-emission" rate for an ac- 
creting white dwarf. The dot-dashed line - accretion rate. 
The dashed line - MEdd, first for a white dwarf, then for 
a neutron star. At point A the He mantle of the donor 
is exhausted, at point B the white dwarf collapses into a 
neutron star. 



In 



contradict the measured Ne/O ratio in 4U 1626-67 ( |5chulz 



model for Lhc population of compact stars in 



et al. 2001): 0.22 ±0.15. 



the Galaxy (Nelemans et al. 2001) the birthrate of ONe 
white dwarfs with M > 1.2 Mq accompanied by "hybrid" 
white dwarfs with 0.35 < M < 0.5 M Q which get into 
contact within 10 Gyr is 1.25 x 10 _6 yr — x . This rate is 
within the limits for occurrence of A IC's in the Galax y 
set by nucleosynthesis considerations ( Fryer et al. 1999 ). 

The initial binary may contain e. g. a 0.38 M Q "hy- 
brid" whi te dwarf [for which we know the internal struc- 
ture from [Ibcn fe Tutukov (1985| )1 and a 1.2 Ma oxygen- 
neon one (Fig. ^J). At the be ginning of the RLOF P or b = 
2min. The outer 0.12 M© of the donor consist of He. 
The maximum rate of stable He-burning at the surface 
of a massive white dwarf is about 10 -5 M Q yr -1 (Kato fc 



Hachis i 1999] ) , lower than the mass loss rate by the donor 



in the initial stages of mass transfer. Mass exchange in this 
case may be also stabilised by isotropical re-emission of op- 
tically thick wind from accretor, generated by He-burning 
(Hachisu et al. 1996). Therefore stable helium burning is 
likely to set in on the surface of the 1.2 Mq accreting white 
dwarf. 



The collapse interrupts mass transfer for ~ 1.56 x 
10 8 yr due to the loss of bind ing energy of the dwarf (see 
for details Verbunt et al. 199C| ) . This time span may not be 
sufficient for the decay of the magnetic field of the pulsar. 
When after resumption of the contact the period of the 
system has increased to 41.4 min, the mass of the donor is 
O.O1O1M . 

Until this moment the neutron star in our model sys- 
tem has lived after the AIC for about 350 Myr. If mag- 
netic fields do not decay on long time scale its field will not 
have decayed sufficiently. The accretion of 0.003 Af Q onto 
a neutron star is most probably insufficient for destruction 
of its magnetic field (i f the frozen field and incompressible 
fluid approximations ( |Cheng fc Zhang 199S ) are assumed) . 
Assuming solid body rotation and conservation of angular 
momentum, one gets ~ 1.3 ms for the initial spin period of 
the neutron star. The present spin period of the neutron 
star in 4U 1626-67 is 7.7 s, clos e to the "death-line " of ra- 
dio pulsars for B = 3 x 10 12 G ( [Bhattacharya 1996| ). Thus, 



Like Verbunt et al. (199C), we assume that the white 
dwarf collapses into a 1.26 M Q neutron star after the in- 
crease in mass to 1.44 Mq. This happens after 240Myrs of 
accretion, when P rb~ 2000 s. The mass of the donor is at 



that moment 0.0132 Mq. The Nc/O ratio in the core for 
azeotropic abundance of Ne (0.1 to 0.3, see 2.2) doesn't 



it is conceivable that 4U 1626-67 harbours a neutron star 
which ceased to be a radio pulsar, but didn't spin down 
sufficiently to experience recycling. 

At the stage when our model system resembles 4U 
1626-67, the mass exchange rate is ~4 x 10 -12 M© yr -1 . 
As already noticed before ( Chakrabarty 1998; [gchulz ct"al 



2001), all model M for 4U 1626-67 based on M - R rela- 



tions are considerably lower than the estimates of M based 
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on the observed spin-up rate and simple angular momen- 
tum conservation considerations. It is plausible that the 
secular M is consistent with model expectations, while 
mass transfer rate M inferred from the observed spin-up 
of the pulsar (if correct) reflects the accretion rate from an 
unstable accretion disk ( |5chulz et al. 2001 ). The latter ex- 
periences thermal ionization instability if Af a ^7.4xl0 -10 
and ^S2.5xlO _9 Moyr~ 1 , for a pure C or O disk, respec- 



tively (Menou et al. 2002) 



4. Discussion 

White dwarfs may be considere d as a ternary C/ O/Ne 
ionic mixture. For such mixtures Ogata et al. (1995| ) pre- 
dict the formation of an almost pure Ne core upon so- 
lidification for any C:0 ratio. If this core ultimately con- 
tains all 22 Ne generated in the white dwarf its mass is 
only ~ 0.004 Mq, much less than given by Eq. (]l|) and this 
would require much larger P rb for X-ray binaries showing 
Ne enrichment. 

On the other hand, according to |5egretain (1996 ), if 
a 0.6 Mq white dwarf is considered as a ternary mixture 
with X c = X Q = 0.495, A 2 2 Nc = 0.01, almost all Ne con- 
centrates in a thin layer around M r /M w< \ w 0.7. Then 
spilling of Ne-enriched matter over the neutron star would 
require extremely fine tuning of the model. 

Another intriguing problem concerns X-ray bursts. 
They were reported for 4U 0614+091, 2S 0918-549, and 
4U 1850-08 57. None of the b ursts have shown atypical 
behaviour (Juctt et al. 2001). Carbon flashes at M < 



0.1A/Edd were never studied. Extrapolation of the calcu- 



lations of Cumming & Bildsten (2001) over 2 orders of 
magnitude may be not very relevant. However, it shows 
that for Xq ~ 0.2, characteristic for "hybrid" dwarfs and 
M ~ 10~ n Mq yr _1 , accreted carbon possibly burns sta- 
bly. Thus, the origin of the bursts in systems with any 
low-mass C-rich donors deserves further attention. 



5. Conclusion 

We have shown that within the framework of the white 



dwarfs cooling model developed by Isern et al. (1991 



Bravo et al. (1992 ), [gegretain et al. (1994 ), Hernanz et al 



(1994 ) it is possible to explain the formation of the Ne- 



enriched donors in X-ray binaries if their progenitors were 
low-mass CO white dwarfs with thick He mantles. A neces- 
sary condition is that the white dwarfs had enough time to 
be substantially chemically fractionated. This means that 
the mass transfer had to start several Gyr after the forma- 
tion of the white dwarfs. We estimate that the masses of 
Ne-rich donors in these systems are ^ 0.06 Mq and their 
orbital periods have to be lOmin. The expected Ne/O 
ratios are ~ (0.1 — 0.4). 

Finally we conclude that the X-ray pulsar 4U 1626-67 
could follow an evolutionary path similar to the other Ne- 
rich systems, with the only difference that it used to have 
a massive ONe white dwarf component that collapsed into 



a neutron star relatively recently. We estimate the current 
mass of the donor in this system as 0.01 Mq. 
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